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We discuss the appearance of magnetic and charge instabilities, named respectively metamag- 
■ netism (MM) and phase separation (PS), in systems which can be described by a perturbed Luttinger 

\ liquid. We argue that such instabilities can be associated with the vanishing of the effective Fermi 

£SJ ■ velocity v, which in some cases coincides with a divergence of the effective Luttinger parameter 

K. We analyze in particular an XXZ chain with next-nearest-neighbor interactions in different 
limits where MM shows up and an extended Hubbard model where in turn, PS occurs. Qualitative 
agreement with previous studies is found. 
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I. INTRODUCTION 



. The study of instabilities in low dimensional strongly correlated electron systems has received much attention in 
the last few years. One of the main reasons is that a charge instability phenomenon (phase separation (PS)) often 
, t \ shows up in the vicinity of the superconducting transition in cuprates. In the case of double exchange models for 
manganese oxides that exhibit the "Colossal" magnetoresistance effect Jjj, this charge instability arises close to the 
transition to ferromagnetism and, interestingly, finite size studies of both the two-dimensional realistic model and its 
one-dimensional version display similar features in this respect. The case of magnetic instabilities (metamagnetism 
(MM)) has also received recent attention in connection to the one dimensional antiferromagnetic (AF) XXZ spin 
chain with next-nearest-neighbor interactions (NNN), where it was found that MM arises in a finite region of the 
q ■ phase space ||. 

— i \ Generically, at the point where charge or magnetic instabilities occur, a divergent compressibility or magnetic 
susceptibility arises. This divergence is in turn associated respectively with the coexistence of two phases with 
^ ' different hole concentrations or magnetizations. 

The aim of this paper is to discuss a general way to determine whether an instability could show up in a given 
Q\ ' one-dimensional model, using Abelian bosonization. Our discussion applies to any one-dimensional model describable 
as an integrable model plus perturbations whose effect is to renormalize the Luttinger liquid (LL) parameters, K and 
v. The effect of irrelevant operators is also discussed. Our discussion could be also relevant in the study of certain 2D 
systems which can realize the so called sliding Luttinger phase ||, since a divergence in the ID susceptibility leads 
to a singularity in the 2D one Q . Generically, our approach provides a quick tool to study the tendency of different 
perturbations to produce instabilities in the above mentioned systems. 

As a sample case for magnetic systems, we analyze the XXZ chain with NNN exchange, treating first this last 
interaction perturbatively within a bosonization approach and find qualitative agreement with the results obtained in 
H in finite systems (This problem was also studied in (|]). We also study the opposite coupling limit, which we call 
hereafter "zig-zag" limit, where the system can be reinterpreted as a two-leg zig-zag ladder. In this case, we determine 
£h . a region in the parameter space where MM occurs. 

To analyze the case of charge instabilities we consider the charge sector of the Hubbard model at incommensurate 
fillings, perturbed by a nearest-neighbor density-density interaction V. In this case we find that for small Coulomb 
repulsion U there is a region where the system phase separates which corresponds to negative values of V, in agreement 
with previous studies ||. 



II. GENERAL DISCUSSION 



We consider a generic situation in which the low energy degrees of freedom (which could correspond to either 
charge or spin variables) are described by a Tomonaga-Luttinger liquid in the unperturbed case. The corresponding 
Hamiltonian once interactions are taken into account will generically be of the form (modulo irrelevant terms) 

H + a(d x ^) 2 + f3(dj) 2 + \(d x <f>) , (1) 
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where Hq corresponds to the Tomonaga-Luttinger Hamiltonian 



H = ljdx (vK(dJ) 2 + £(0«,$ a ) • (2) 

The first two terms can be readily absorbed into a redefinition of the LL parameters K — > K e f f and v — > v e f f , 
while the third one changes the chemical potential in the case of charge variables and the magnetic field in the spin 
case. 

Under the above mentioned conditions the compressibility for the charge modes, described by an effective LL with 
parameters K c and v c can be shown to be given by |t]] 

K c 

K OC — , (3) 

V c 

and similarly, the magnetic susceptibility for the spin modes described by an effective LL with parameters K s and v s 
is given by 

X«— • (4) 

V s 

One readily observes that a divergence in these quantities arises either when K7 cs *. or W( CjS ) vanish, ft should be 
pointed out that these two things could happen simultaneously, but this is not always the case and hence PS or 
MM instabilities are to be identified with the vanishing of the effective velocity. We argue that this is true, provided 
that the quantities Kv and K/v remain positive definite and irrelevant perturbations do not change the large scale 
behavior. In Ref. JsJ an attempt to characterize these instabilities in one dimensional systems exhibiting a universal 
character described by a Tomonaga-Luttinger model was made. In this paper the authors identified the divergence of 
the above mentioned thermodynamic quantities with a divergence in the so called Luttinger parameter K. Although 
this interpretation led to a consistent analysis for the cases studied in we argue that a more general criterion 
consists in identifying the instability regions with those where the velocity vanishes. This last statement can be 
understood as follows: by analogy with the unperturbed XXZ chain in a magnetic field, we see that when the Fermi 
velocity goes to zero we approach a ground state of FM nature (and K does not necessarily diverge, though this 
happens for the particular case of the XX Z chain for A = — I and zero field ^|). Then the magnetic susceptibility 
diverges as I/w e // and, if this happens before reaching saturation, the magnetization curve as a function of the applied 
magnetic field presents a jump. Once we have shown that there is a jump in the magnetization curve, for some critical 
value of the applied magnetic field, h c , we can conclude that two different values of the magnetization will coexist at 
this point 0. 

In the next two Sections we study two sample cases where this general discussion applies, the XX Z chain with NNN 
interactions where MM has been shown to occur Jp|, |i| and the extended Hubbard model with nearest-neighbour 
interactions V, where PS appears for negative V 



III. THE XXZ AF CHAIN WITH NNN INTERACTIONS 



The Hamiltonian of the XXZ AF with NNN interactions in a magnetic field is given by 



N 

Hxxz — J {SfSf+i + SfSf +1 + AS*S* +1 ) + 

i=l 

N N 

J' ]T (SfSf +a + SfSf +2 + AS?S? +2 ) -h^S?. (5) 

i=l i=l 

The large scale behaviour of the XXZ chain can be described by a U(l) free boson theory with Hamiltonian (g). 
The field 0* and its dual (jf are given by the sum and difference of the light-cone components, respectively. The 
constant K governs the conformal dimensions of the bosonic vertex operators and can be obtained exactly from the 
Bethe Ansatz solution of the XXZ chain (see e.g. |[o| for a detailed summary). We have K = 1 for the SU(2) 
symmetric case (A = I) and is related to the radius R of (To) by K^ 1 = 2nR 2 . In (g) v corresponds to the Fermi 
velocity of the fundamental excitations of the system. 

In terms of these fields, the spin operators read 
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S x ~ -}=d x (f> + a : cos(2k F x + V2^(f>) : , (6) 

v 27r 2 



5' 



± 



(_1)« : gWEF* (bcos(2k F x + VW>) + c) : , (7) 



where the colons denote normal ordering with respect to the groundstate with magnetization (M). The Fermi 
momentum hp is related to the magnetization of the chain as k F = (1 — (M))w/2. The effect of an XX Z anisotropy 
and/or the external magnetic field is then to modify the scaling dimensions of the physical fields through K and 
the commensurability properties of the spin operators, as can be seen from (^), (Q). The constants a, b and c were 
numerically computed in the case of zero magnetic field [jllj (see also jl2| ) . 

In what follows we study both the weak coupling (J'/ J <C I) and the zig-zag (J'/ J ^> 1) limits. 

i) J'/ J < 1 limit 

Using (^|) and (Q) the NNN interaction term in the bosonic language reads 

H N nn = a J dx (gi (d x (f)) 2 + g 2 (d x <f?J + A d x <fi + A' cos(V^4>)j . (8) 
where a = J' / J and <?i,2, A and A' depend on A and the non-universal constants a, b and c as 

gi = A ( 2™ 2 cos(4Mj + 27r6 2 cos(4fc F ), (9) 
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//,>::■ -1,7 ( ^ - ^-COSU/.", ) ) . ( 10) 



A=— LoA(M), X' = abc. (II) 



and 

A = - 

2tt 

The first two terms in (||) have the effect of renormalizing both the compactification radius and the Fermi velocity 
in the following way 

v 2 eff = (v/K + 2a 9l ) (vK + 2ag 2 ) , K 2 eff = -——2L , (12) 

Now we can make contact with the analysis of Ref . ||] : the MM region is identified within our approach as the set 
of phase space points in which the susceptibility diverges for < M < M sat . 

One immediately observes from ( |l2] ) that, for zero magnetic field and without taking into account the effect of 
the A' perturbation, both v e ff and K~^f vanish simultaneously and that this happens when (v/K + 2agi) = (see 
Fig. 1). However, this situation can change for two reasons: first, an extra term arises from the NNN perturbation, 
which rcnormalizes the external magnetic field h and hence changes both v e ff and K e ff through the Bethe Ansatz 
equations, generically in a way which removes the above mentioned simultaneity. Second, if one goes beyond the zero 
loop order, both K e ff and v e ff will renormalize in a different way due to the A' term p3| . We will discuss this issue 
again in the context of charge instabilities in the next Section. In particular, the boundary between the MM and 
the FM phases is obtained as the set of points in which v e ff and/or vanish/es for h = 0, and that between 

the MM and the AF phases as the set of points in which v e f / ^ and ^ for all values of the magnetization 

(M) < M sat . 

Let us focus on the MM-FM boundary, which is easily obtained using (|) and ( 12|) . In this case, we can use the 
numerical values for the non- universal constants a, b and c appearing in (|) and (f?) obtained in [p"l| . The boundary 
obtained in this way agrees qualitatively with that obtained in Q (see Fig. 1). The lack of quantitative agreement 
is presumably due to the perturbative treatment of the NNN interactions, which could be improved by considering 
higher loop contributions in a renormalization group analysis from the bosonization side and due to finite size effects 
from the numerical one. However, our main aim is to discuss the appearance of instabilities in a generic and simple 
way. 



3 



One should be careful also about the regime of validity of this approach, since due to the renormalization of 
the Luttinger parameter K e ff, the scaling dimensions of the many discarded irrelevant perturbations change, and 
nothing prevents one of these to become relevant. In fact, by analyzing the scaling dimension of the leading irrelevant 
A' perturbation, we observe that it reaches the limiting value 2 at a c = ^Kfa-gl) ' an< ^ ncnce our approach ceases 
to be valid for a > a c . This critical line separates the massless regime from a massive one, and our results compare 
qualitatively well with those obtained in (dashed line in Fig. 1). The study of the AF-MM boundary is in this 
case more involved due to the appearance of the non-universal constants a, b and c in the bosonized operators which 
are not available for non-zero magnetic field. (The estimation of the field dependence of these constants for < A < 1 
but its evaluation is more delicate on the ferromagnetic side A < [pi.) 
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FIG. 1. Ground state phase diagram for the XX Z chain with next-nearest-neighbor interactions in the limit J'/J <C 1 
showing the transition lines between the Luttinger liquid regime and the ferromagnetic one. Square dots are reproduced from 
[2] and the solid black line corresponds to our theoretical prediction. The dashed line indicates where the leading irrelevant 
operator becomes marginal. 



ii) Zig-zag limit ( J'/J > 1) 

In 4- l-i n o lmd nm r/lmf ArrvvrtI 4- 



In this limit we reinterpret the Hamiltonian (j^) as a two leg zig-zag ladder in the weak interchain coupling limit 
fl7| p2| . In this description one represents each of the chains by one free compactified £7(1) boson. Thus, the whole 
ladder is represented by two bosons <f>\, <j>2, each governed by an action given by (||), plus the perturbative terms 
arising from the interchain zig-zag coupling. 

This system, in the presence of an external magnetic field, has been studied in |22| , where it was shown that 
for non-zero magnetization a relevant interaction gives a mass to the "relative" field <j) re i = (</>i — </>2)/2, while the 
"diagonal" field, 4>diag = (</>i + <fe)/2, remains massless. 
The effective Hamiltonian governing the large scale behavior of ifidiag is then given by 

H ef f = ~ J dx {v e ffK eff {d x 4>diag)' 2 + J^JJ {9 X 4>diag) 2 + \d x (j)diag^j , (13) 



where, calling (3 = l/a, A = /3A(Af)/V27r and v e ff and K e ff are given by 



■l /=w 2 ll + ^l), Kl // =^ll + ^l) . (14) 



PAK\ 2 2 / (3AK \ 1 
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There is however an extra term that mixes the fields <fi re i and 4>diag 



= -A' / dxd x (/)diag sin(V^TT(f>rel) 



(15) 



The effect of this term was studied in Ref . |20j where it was shown that it gives rise to a spin nematic phase close to 
A = 0. The main point for our analysis is that the diagonal field is still described by a LL, and we expect the same 
picture to apply for —1 < A < 0. For this reason, we do not take into account this last term in what follows. 

One then observes from ( fbif ) that the non-universal constants do not appear in the perturbing terms, and hence we 
can obtain in this case both the FM-MM and MM-AF boundaries. The results are presented in Fig. 2. 
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FIG. 2. Qualitative ground state phase diagram for the XXZ chain with next-nearest-neighbor interactions in the limit 
J'/J»l. 



IV. HUBBARD MODEL WITH NEAREST NEIGHBOR INTERACTIONS 



We consider now the Hubbard model, defined as 



(16) 



where 



(17) 



In the absence of an external magnetic field it has been shown that this model presents charge-spin separation. 
It has also been shown (c.f. ^3|) that the large scale behavior of the spin and charge degrees of freedom can be 
described by two decoupled boson field theories with dynamics governed by the Tomonaga-Luttinger Hamiltonian. 
The parameters K and v can in each case be exactly obtained for all values of /i and U via numerically solving the 
Bethe Ansatz equations in [^i) . Approximate expressions for the velocity of the charge sector in the small and large 
U regimes are given in ]25( ] . 

The addition of a density-density interaction between nearest neighbors leads to the so called Extended Hubbard 
model, whose Hamiltonian is given by ( |l6| ) plus the term 



5H = VJ2>, 



IjTlj+l, 



(18) 



where 



(19) 
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It is a simple matter to show that, modulo irrelevant operators, the effect of the perturbation (|T^) is to renormalize 
the parameters [i and U as follows: 



fi + V, 



U^U + 2V. 



(20) 



Therefore, the low energy behavior of the charge sector of the Extended Hubbard model is that of a Luttinger liquid 
with parameters 



K eff =K(U + 2V,n + V), 



"eff 



v(U + 2V,fi + V), 



(21) 



where K{U,y) y v(U, /i), are given by the exact Bethe Ansatz solution of the non-Extended Hubbard model (|l^). 
This allows for the determination of the divergences in the compressibility k, as given by (|^), which we identify with 
the zeros of the effective velocity. 

We have only considered small values of V, since our approach is a perturbative one. Within this restriction, we 
have found no singularities, and therefore no PS, in the large U regime. On the other hand, regarding the small U ,V 
region, in which the effective velocity is given by 



V eff 



v(U + 2V, 11 + V) ~ ! (l + v 7 ! + 2(m + V)) (a* + V) + 



U + 2V 
2tt 



(22) 



we have found that the roots of the above expression present the same qualitative behavior as that described in [[| . 
It should be stressed that K e f j remains finite in the region where ( ]22] ) vanishes. 



} 

.75 
0.5 
0.25 










.5 -0.25 

-0.25 

-0.5 

-0 .75 

-1 


0.25 


. 5 


.75 


1 1 . 



FIG. 3. Phase separation lines in the extended Hubbard model for V = 0, —0.3, —0.6, —0.9 from left to right. 
Following this approach, one could also study the appearance of instabilities in the magnetic sector. 

V. CONCLUSIONS AND DISCUSSION 

We have discussed the appearance of MM and PS in systems that can be described as perturbed Luttinger liquids. 
Specifically, we have studied the XXZ spin chain with NNN interactions and the extended Hubbard model with 
nearest-neighbor density-density interactions. We have found FM, AF and MM regions in the zig-zag limit of the 
XXZ-NNN spin chain. In the weak coupling limit, qualitative agreement with previous results was found for the 
MM-FM transition. We were not able to complete our search for MM due to lack of specific numerical data, an issue 
to be discussed elsewhere. The instabilities were in all cases identified as the roots of the inverse susceptibility (|4|). 

Concerning the extended Hubbard model, according to our analysis, it does not present PS in the large U and small 
V limit. On the other hand, in the small U, V limit instabilities do show up, through the roots of the effective Fermi 
velocity, while the effective Luttinger parameter remains finite. 
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We argue that all the instabilities studied are associated in general with a vanishing effective velocity V e ff. As a 
matter of fact, it is the roots of v e ff that give a divergent compressibility of the extended Hubbard model, though in 
the other cases studied, this coincides with a divergence of the effective Luttinger parameter K e f / . This should not 
be regarded as the actual hallmark of the instabilities: although our treatment of the XXZ-NNN spin chain gives at 
the same time v e ff — — 0, this is only true in the absence of an external magnetic field and moreover ceases to 

be valid beyond the zero loop order since both parameter renormalize in a different manner [ fL3| . It should be stressed 
that our approach gives a simple way to qualitatively analyze instabilities in generic charge and magnetic systems 
provided they can be described as perturbed Luttinger Liquids. Henceforth, it provides a quick tool to study whether 
different perturbations could produce such instabilities. 
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